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Abstract. Solitary wave structures observed by the Ulysses spacecraft 
in the solar wind were analyzed using both inverse scattering theory as well 
as direct numerical integration of the derivative nonlinear Schrodinger (DNLS) 
equation. Several of these structures were found to be consistent with soli- 
ton solutions of the DNLS equation. Such solitary structures have been com¬ 
monly observed in the space plasma environment and may, in fact, be long- 
lived solitons. While the generation of these solitons may be due to an in¬ 
stability mechanism, e.g., the mirror instability, they may be observable far 
from the source region due to their coherent nature. 
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1. Introduction 

For several decades, structures in the interplanetary magnetic field, called “holes,” “an¬ 
tiholes,” “bumps,” and “dips,” have been observed. One of the first published observations 
was of data taken by the magnetometer on the Explorer 43 spacecraft (also known as Imp 
I), where Turner et al. [1977] saw localized regions of very low intensity magnetic field 
(less than 1 nT). They, along with Burlaga and Lemaire [1978] tendered a theoretical 
explanation: these regions were “stationary nonpropagating equilibrium structures.” To 
deduce this, they used a kinetic boundary layer theory [Lemaire and Burlaga, 1976] to re¬ 
produce the magnetic field profiles that were observed. Another, more recent, explanation 
for these structures is that they are freak waves (or rogue waves). Laveder et al. [2011] 
has shown that rogue waves appear in randomly driven, dissipative DNLS systems, but 
it is not clear how to include these effects in an analysis of the solar wind. In addition, 
Ruderman [2010] has shown that freak waves can arise in the DNLS context due to both 
linear and nonlinear dispersive focusing. 

In subsequent observations, it has become accepted that, because these structures are 
spatially near regions that are unstable (or at least near marginal stability) against the 
mirror instability, they are mirror mode structures. As Winterhalter et al. [1994] put it, 
“there are indications that the holes may have been remnants of mirror-mode structures 
created upstream of the points of observation.” In the two years of Ulysses observations 
(1991-1992), Winterhalter et al. [1994] found about 30 events per month in which the 
magnetic field did not change direction across the hole, a type that they called “linear 
holes.” Hamilton et al. [2009] has made the case that these types of holes are well-suited 
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to analysis using the derivative nonlinear Schrodinger (DNLS) equation and the accom¬ 
panying inverse scattering transformation (1ST). In addition, using Wind magnetometer 
observations, Stevens and Kasper [2007] also concluded that “many magnetic holes are 
remnants of the mirror-mode instability.” 

Rees et al. [2006] catalogued about 30 events in the Ulysses data set, i.e., localized 
increases in the magnetic held strength that had “the characteristics of a soliton.” They 
then compared them with the bright solitons of Baumgdrtel [1999] which were solutions 
to the Hall MHD equations, and were similar to the so-called “one-parameter” soliton 
solutions to the DNLS equation. However, while the properties of the comparison was 
suggestive, the agreement between the observations and the analytic solution was not 
satisfactory. 

We will show below that the problem was twofold. First, for the event that they 
analyzed in detail (shown in their Fig. 8), an 1ST analysis using the transform for the 
DNLS equation revealed no eigenvalues, and hence no solitons. And second, even if the 
overall prohle shows an increase in the magnetic held strength, this does not preclude the 
possibility that the primary soliton embedded in the prohle is a dark soliton. 

In this paper, we re-analyze several of the magnetic structures observed by Ulysses and 
catalogued by Rees et al. [2006], but we use the suggestion of Hada et al. [1993]. That 
is, we attack each prohle from the perspective of inverse scattering theory and the DNLS 
equation. Each magnetic structure is analyzed in two ways. First, the 1ST technique is 
used to obtain an analytic soliton solution that is embedded within the observed magnetic 
held. This method searches for eigenvalues of the associated linear eigenvalue problem, 
each of which corresponds to either a one-parameter or a two-parameter soliton. Second, a 
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direct numerical solution to the DNLS equation is obtained using the observed magnetic 
held structure as an initial condition. A comparison is then made between these two 
methods. It is found that in the long-time limit, after any dispersive component has 
radiated away, the numerical solution matches the analytic soliton well. This result implies 
that DNLS solitons are the main component of these magnetic structures. 

Of course, it would be more accurate to numerically solve the MHD equations or the 
two-huid equations rather than the DNLS (to evaluate the time evolution of the initial 
magnetic held prohle). However, since Ulysses captures only a snapshot of the magnetic 
held, it is more enlightening to compare with the DNLS time evolution, given that an 
eigenvalue (or eigenvalues) was found via the 1ST, signifying the presence of a soliton. 


2. The derivative nonlinear Schrodinger (DNLS) equation 

The DNLS equation is a general nonlinear wave equation that describes solitons in 
plasmas under widely varying conditions. While the standard derivation is for describing 
quasi-parallel Alfven waves, in various limits it reduces to the Korteweg-de Vries (KdV) 
equation and the modihed KdV equation [Kennel et ai, 1988]. These can be used to 
describe ion-acoustic solitons, for example [Washimi and Taniuti, 1966]. 

A useful form of the DNLS equation is given by Baumgdrtel [1999], 


db d 
dt ^ dx 




1 d‘^b 

2 dx‘^ 


0 , 


( 1 ) 


where b{x, t) = by + ibz and 

4cos20-/3' 

In this case, the wave propagates in the x direction, where x is scaled to the ion inertial 
length c/ujpi, and Eq. (1) is written in a frame moving in the x direction with the 
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intermediate Alfven speed ca cos 6 relative to the background plasma. In addition, the 
magnetic held b has been scaled to the background magnetic held strength Bq away from 
the soliton. Also, Upi is the ion plasma frequency, c is the speed of light, (3 = cl/c\ is the 
plasma “beta,” Cg is the plasma sound speed, ca = Bq /is the Alfven speed, and p 
is the plasma mass density. With a suitable redehnition of x and t it is possible to set 
the coefficient a equal to unity (see Baumgdrtel [1999] and Mjplhus and Hada [1997]), 
and we use this redehnition to simplify the numerical calculations. Other versions of the 
DNLS equation can be found in Kennel et al. [1988] and Hamilton et al. [2009]. The fact 
that bo is nonzero rehects the fact that the soliton is not propagating exactly parallel to 
the background magnetic held, and therefore the boundary conditions are “nonvanishing” 
[Kawata and Inoue, 1978]. 

Several analytic solutions to this equation have been found, and the most well known fall 
into two classes, the so-called one-parameter (either dark or bright) and two-parameter 
solitons. (For explicit expressions for these solutions, see, for example, Mjplhus and Hada 
[1997] or Baumgdrtel [1999] or Hamilton et al. [2009].) These single soliton solutions 
can also be found using the 1ST technique, which uses the fact that some nonlinear 
wave equations, e.g., Eq. (1) with solution b{x), are associated with Sturm-Liouville-like 
problems where b{x) takes on the role of the potential function. For a given b{x), the 
scattering data results in eigenvalues that correspond to soliton solutions of the nonlinear 
wave equation. For the DNLS equation, the 1ST is a 2 x 2 eigenvalue problem, and the 
entire procedure is known as the “AKNS scheme” [Ablowitz and Segur, 1981]. As shown 
by Kawata and Inoue [1978], this eigenvalue problem can be written in the form 

$. = 7^ ■ $, (3) 
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where 


A) = 


011 012 
021 022 


and 


D = X 


—i\ b{x,t) 
b{x,t) iX 


( 4 ) 


( 5 ) 


and where b{x, t) is the wave prohle of the magnetic held and A is the eigenvalue. The one- 
parameter solitons mentioned above are characterized by real A, and the two-parameter 
solitons have complex A (the real and imaginary parts of A are the two parameters). A 
general magnetic held prohle, such as the Ulysses prohles examined here, will be composed 
of solitons (denoted by discrete eigenvalues) as well as dispersive waves that are part of 
the continuous spectrum. 

To search for eigenvalues of a particular prohle b{x), the relation between the eigen¬ 
functions is conventionally given through a scattering matrix 


$ {x, t] A) = $+(a;, t; A) • J , (6) 

where the elements Sjj(A) are functions of the eigenvalue A, and are the eigenfunctions 
of Eq. (3) dehned through their asymptotics as discussed in the Appendix. Given b{x), 
and using the dehning asymptotics of as boundary conditions, a numerical solution 
can be found for $ which can then be used to obtain Sn. As is standard in scattering 
matrix theory, the condition Sii(A) = 0 determines the set of eigenvalues {A„} for which 
the eigenfunctions are bounded both to the left and right. The existence of a discrete 
eigenvalue is essentially the dehnition that a soliton is embedded in b{x) [Drazin and 
Johnson, 1989]. In addition, dispersive wave solutions to the DNLS (corresponding to the 
continuous spectrum of eigenvalues) can also be present in b{x), or they can be the only 
component. 
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The numerical technique that we use to solve the DNLS equation is a two-step Lax- 
Wendroff method for the nonlinear term with a backward-time, centered-space method 
for the diffusive term. The boundary conditions are periodic, and a large domain is used 
to ensure that dispersive waves do not have time to propagate around the domain and 
interfere with the soliton. 

3. Ulysses events 

Several detailed studies of Ulysses observations of magnetic holes have been made [Tsu- 
rutani et ai, 1992; Winterhalter et ai, 1994; Franz et ai, 2000; Rees et al, 2006]. Because 
of the large speed that the spacecraft is moving through the solar wind plasma, and the 
small Alfven speed at which the structures are moving through the plasma, each obser¬ 
vation is effectively a snapshot of the magnetic held at one instant of time. 

The vector helium magnetometer on Ulysses provided one-second averaged magnetic 
held measurements, reported in RTN (Radial-Tangential-Normal) coordinates [Balogh 
et al, 1992]. An example of the raw data, from 17 Jul 2002, is shown in Fig. 1(a). 
The characteristics of the plasma environment during the magnetic held observation, pro¬ 
ton temperature, plasma velocity, and density, were obtained by the SWOOPS (Solar 
Wind Observations Over the Poles of the Sun) experiment [Bame et ai, 1992]. Typical 
values were (3 ~ 0.1, plasma speed Vpiasma ~ 500 km/s, and ion inertial length c/upi ~ 10^ 
km. 

Rees et al. [2006] surveyed the entire Ulysses database looking for soliton-like structures 
in the magnetic held. Their criteria were “short duration increases in the magnetic held 
magnitude and an associated rotation in the magnetic held direction.” However, since 
the magnetic held vector was in the same direction on both sides of the structure, these 
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are “linear” in the terminology of Winterhalter et al. [1994], Rees et al. [2006] found 33 
events, and while they appeared to be soliton-like, the quantitative comparison with the 
analytic form of a bright, one-parameter soliton was not “adequate.” Here, we re-examine 
some of the events found by Rees et al. [2006] and hnd that the 1ST results in eigenvalues 
for many of the events. 

The events we have chosen for careful study are shown in Table 1. The hrst event listed 
in Table 1, on 21 Feb 2001, was studied in detail by Rees et al. [2006] (it is shown in 
their Figure 8). Because it consisted of a magnetic compression, they attempted to model 
it with a bright soliton solution of the DNLS equation [Baumgdrtel, 1999]. However, as 
can be seen from Fig. 8 of Rees et al. [2006], the £t was not very good. Our explana¬ 
tion for this poor agreement is that there is no DNLS soliton embedded in the profile, 
which is verihed by the lack of a discrete eigenvalue found by our scattering analysis. 
As mentioned above, a discrete eigenvalue determines the presence of a soliton, and the 
lack of any eigenvalue implies that there are only dispersive wave solutions to the DNLS 
in the magnetic held profile. The subsequent evolution of that magnetic held profile, as 
determined by a numerical solution of the DNLS, shows that the initial pulse collapses 
into a wave train. Another reason for the lack of a DNLS soliton is that this wave was 
propagating almost perpendicular to the background magnetic held. This is a regime 
where the DNLS equation is thought not to be applicable for low P plasmas. See, for 
example, the discussions in Baumgdrtel [1999], Buti et al. [2001], and Ruderman [2002]. 

The second event listed in Table 1, on 17 July 2002, is shown in Figs. 1-3. Figure 1 
(a) shows the raw magnetometer data, while 1 (b) shows the rotated, smoothed, and 
scaled held components. A minimum variance analysis [Sonnerup and Scheible, 2000] was 
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performed to determine the orientation of the structure, and the x-direction was scaled 
in order to use Eq. (1) with a = 1. Finally, the prohle was smoothed using exponential 
smoothing in order to simplify the root-hnding algorithm that searched for eigenvalues. 

Figure 2 shows the magnetic held structure from Fig. 1(b), which is the initial condition, 
along with a large-amplitude one-parameter dark soliton with eigenvalue Ai = 0.223. 
There was a second real eigenvalue found by the scattering analysis, but this had a very 
small amplitude since its eigenvalue was close to unity, A 2 = 0.994, and in addition, there 
were dispersive waves. That the Ai = 0.223 soliton is the most signihcant can be seen 
from the energy in the magnetic held 

/ OO 

{\b\‘^-bl)dx. (7) 

-OO 

This integral can be analytically evaluated for solitons [Sanchez-Arriaga, 2010] and pre¬ 
dicts that the total energy of the two dark solitons is Ei = —5.38 and E 2 = —0.438, 
respectively, and this was conhrmed numerically. In addition, the energy of the dispersive 
waves, which can only be found numerically, is Edisp = 3.28. This shows that the second 
soliton has minimal energy compared with the dispersive waves, and once they propagate 
away the solitary structure is ehectively a single soliton. 

This can be seen in Fig. 3, which shows the observed magnetic held and the DNLS 
soliton at time 55.2 later, where Qp is the proton cyclotron frequency. The Ulysses 
magnetic held has been propagated forward in time with a numerical solution of the DNLS, 
and compared with the analytic time-dependent dark soliton solution. Note that the 
soliton matches the observed magnetic held much more closely at this later time, because 
the dispersive wave component has propagated away toward the edges of the simulated 
domain. Also, this analysis answers the question of why a dark soliton represents a 
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magnetic increase so well. The reason is that (in this case) the compressive magnetic held 
prohle consists of a large-amplitude dark soliton, a small-amplitude dark soliton, and a 
dispersive wave component with signihcant energy, as shown above. The stability of dark 
solitons allow them to remain embedded in the magnetic held, while any bright solitons 
decay away. 

4. Discussion 

It is instructive to note that all of the events in Table 1 are characterized by multiple 
eigenvalues. This simply means that the wave form of the magnetic held is comprised of 
multiple solitons. Of course, the combination of solitons is not linear, but nonlinear, and 
general methods to obtain so-called “iV-soliton solutions” can be found in Ablowitz and 
Segur [1981] and Drazin and Johnson [1989], for example. In certain cases, such as the one 
examined in detail in Figs. 1-3, only one soliton contributes signihcantly to the magnetic 
held prohle. In general, however, the nonlinear iV-soliton solution must be compared with 
the observed prohle. It is important to note also that the existence of a discrete eigenvalue 
for an arbitrary prohle, b{x,t), does not necessarily mean that the prohle is a soliton. It 
simply means that there is a soliton embedded in the prohle, possibly along with another 
soliton (or solitons) or a dispersive wave component, or both. 

As long as the underlying nonlinear dynamics are captured by the DNLS, then solitons 
are the most fundamental way of representing the wave dynamics. In addition, if the 
plasma conditions are consistent with dissipation or nonlinear Landau damping or even 
random density huctuations, then it is plausible that the multiple solitons have a common 
origin. That is, the damping of a bright one-parameter soliton will inevitably give rise 
to the formation of a two-parameter soliton which will subsequently damp away itself 
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and form a train of dark one-parameter solitons in the process [Sanchez-Arriaga, 2010]. 
This naturally leads to the question of how the solitons that we have identihed, moving 
with different speeds, are able to arrive simultaneously from a distant source location. 
Of course, within the context of the pure DNLS equation, they cannot. However, the 
physical effects just listed are not included in the DNLS, and it is a challenge for future 
work to determine if a combination of these effects could maintain a sufficiently coherent 
wave packet comprised of mostly dark solitons. 

As stated by Hada et al. [1993], this method of analysis is “more efficient than the Fourier 
transform” when analyzing nonlinear waves. They go on to point out that there are several 
other advantages to the inverse scattering analysis, including the fact that both the past 
history and future evolution of the nonlinear waves is contained in the scattering data. 
With a single spacecraft, unfortunately, it is at best difficult to prove these connections. 
However, both the current Cluster H set of satellites and the future Magnetospheric 
Multiscale mission provide opportunities for such analysis. In fact, Stasiewicz et al. [2003] 
showed that there was good agreement between Cluster observations in the magnetopause 
boundary layer and magnetosonic solitons, and Trines et al. [2007] saw good agreement 
between simulations and Cluster observations at the magnetopause of solitary electrostatic 
structures. However, we have shown here that a deeper understanding of nonlinear waves 
can be obtained using the inverse scattering analysis to obtain eigenvalues, in addition to 
a direct comparison between observations and theory. 

Appendix: Inverse Scattering Transform of the DNLS 

For the case of nonvanishing boundary conditions (i.e., when 6 —)■ 6o as a; —)■ ±cxd), 
Kawata and Inoue [1978] showed that the linear eigenvalue equation in Eq. (3), $ 3 , = 


DRAFT 


February 24, 2015, 4:38am 


DRAFT 



WHEELER ET AL.: ULYSSES ALFVEN SOLITONS 


X - 13 


has the asymptotic eigenfunctions 

t] A) — )■ T(A, ()J{Ax) when x —)■ ± cxd , (Al) 

where 

r(A.o=(A''c V) • 

and 

J(Ax) = ( % ) , (A3) 

and where A = X( and ( = a/A^ — b^. Since ( and A are complex, the root chosen to 
express is that which results in a bound eigenfunction, i.e., Im(A) > 0. Then, as 
is mentioned in the text, the scattering matrix S is dehned by ■ S, and the 

requirement that Sn = 0 results in the set of eigenvalues for which <h is bounded. 
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Table 1. Properties of selected events. Columns are date and time (UT), V is the s/c velocity 
relative to the surrounding plasma, /3 is the plasma beta, 0 is the angle between the background 


magnetic held and the layer normal k, and va is the Alfven speed of the surrounding plasma. The 


last two columns list the number of one and two-parameter eigenvalues found for each prohle. 


Date 

Time 

(UT) 

D 

(10^ m/s) 

/3 

0 

Va 

(km/s) 

one- 

parameter 

two- 

parameter 

2001 Feb 21 

03:35 

5.19 

0.498 

93.1° 

38.99 



2002 Jul 17 

08:15 

4.96 

0.085 

32.8° 

98.28 

2 


1991 Jul 02 

21:10 

3.85 

0.153 

31.3° 

88.86 

1 

1 

1991 Jan 05 

04:34 

2.68 

0.058 

26.6° 

103.8 

3 

1 

1991 Apr 01 

04:57 

5.61 

0.123 

15.1° 

83.76 

5 

3 

1992 Jan 14 

19:56 

5.15 
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Figure 1. Ulysses magnetometer data from 17 July 2002, around 08:15 UT. (a) Field compo¬ 
nents as a function of time expressed in RTN coordinates, (b) Field components as a function 
of the scaled and rotated x, after the following transformations have been applied: i) rotation 
using the minimum variance transformation, ii) scaling to the asympototic value of and hi) 
rotation about the new x direction to eliminate the asymptotic value of bz- The eigenvalues that 
resulted from the minimum variance analysis [Sonnerup and Scheible, 2000] were 33.88, 15.75 
and 1, scaled to the minimum eigenvalue. The minimum variance direction, x, makes an angle 
of 6^ ~ 32.8° with respect to the background magnetic held. 
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Figure 2. Solid lines: components of b from Fig. 1(b). Dashed lines: analytic dark soliton 
solution of the DNLS (eigenvalue of A ~ 0.223) at t = 0. Note that the components of 6 at t = 0 
have been padded and the domain increased so that the dispersive waves that are also present 
in the prohle do not have time to propagate around the domain and interfere with the soliton. 
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Figure 3. Same as Fig. 2, but at t = 55.2 The observed magnetic field (solid lines) 
have been evolved using the numerical solution of the DNLS, while the analytic solution (dashed 
lines) has the time dependence included. The analytic solution travels at a speed between the 
Alfven speed and the magnetosonic speed. The particular time was chosen to allow the dispersive 
waves enough time to propagate away from the soliton, but not enough time to wrap around the 
domain. 
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